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SUMMARY 


This report concerns the problem of constructing solutions for 
transonic flows over symmetric airfoils. The aspect of the problem 
emphasized is, of necessity, not how to form a solution for compressible 
flow but how to simplify the initial phase of the problem, namely, the 
mapping of the incompressible flow. In the case of the symmetric 
Joukowski airfoil without circulation, the mapping is relatively simple, 
but the coefficients in the power series are difficult to evaluate. As 
a result, the problem requires simplification. Instead of the exact 
incompressible flow past the airfoil, an approximate flow is used, which 
is derived from a combination of source and sink. This flow differs 
only slightly from the exact one when the thickness is small. By the , 
same method, the flow with circulation is also considered. 

After the incompressible-flow functions are approximated in this 
fashion, the numerical calculation of the corresponding compressible 
flow, by the holograph theory, does not present any essential difficulty. 

INTRODUCTION 


The problem of transonic potential flows has been considered in two 
previous reports (references 1 and 2). The object was, in the first 
place, to construct a solution for a closed body, and, secondly, to devise 
a method by which the flow can be easily calculated. The method / of con- 
structing a solution is essentially a method of analytic continuation. 

That is,' a complex potential, which is known as a function of the complex 
velocity, is represented by a Taylor series about the origin of the holo- 
graph plane and by another, series in an annular region. Both series are 
convergent in their respective domains and agree on the common circle of 
convergence. In the interior’ region, a solution for compressible fluid 
can immediately be obtained by replacing the proper particular integrals 
in the Taylor series. The new solution does not affect the radius of 
convergence and, as Mach number tends to zero, reduces to that of the 
incompressible flow. The solution for the annular region cannot be 
formed by simply replacing the corresponding particular integrals in 
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the outside series of the incompressible flow, because the solution so 
obtained does not agree with the inside series on the circle of con- 
vergence and, consequently, does not represent the same function. 

The problem is then not so simple as it first appears to be. In 
order to make both series represent one solution, the idea of retain in g 
the coefficients in both series of the incompressible flow must be 
abandoned. Instead, both the coefficients as well as the form of the 
solution are modified to satisfy the conditions on the circle of con- 
vergence (reference 2). A different approach to the same problem is 
that due to Lighthi-11 and Cherry (references 3 and 4). Instead of 
choosing in advance the form of the outside series, they first transform 
the Taylor series into a double series, which, on interchanging the 
order of summation, can be summed to yield a single series. As a result, 
the series is valid anywhere except at the singularities in the domain 
considered, and analytic continuation is automatically accomplished. 

The continued outside series consists of two series; one is derived 
from the outside series of the incompressible flow and the other is a 
Taylor series. 

The status of the hodograph method is then this: If a complex 

potential is given in series form,, a corresponding solution for com- 
pressible fluid can be constructed, the body shape being unknown until 
the end of the calculation. It is , apparent that the fundamental problem 
of determining the flow pattern about a given body remains unsolved. 
Moreover, even the transformation of the complex potential often poses 
many practical difficulties. Nevertheless, the hodograph method is 
still the only one available by which an exact solution may be sought. 

For this reason, there is still ground for further exploration along 
this line. 

It is the object of this report to examine a case of more practical 
interest than that previously considered, in order to discover possible 
means of simplification. As a first attempt, a symmetric Joukowski air- 
foil is studied. In this case, if the rule of transforming incompressible 
to compressible flow is strictly followed, the calculation becomes 
extremely laborious. However, if a special case of small thickness - 
which is also of practical interest - is considered, the functions can 
be simplified to such an extent that numerical calculation is practical. 
The process is, in a sense, approximate. But as far as the compressible 
flow is concerned, the merit of a method lies only in whether it yields 
a good aerodynamic body or not. It -has never been demonstrated that an 
exact incompressible potential function gives a better result. From this 
point of view, this is just another way of choosing a complex potential. 
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Finally, it must be emphasized that the present investigation is 
only the beginning, as far as practical calculations for useful airfoil 
profiles are concerned. . Taking the approximate function as a guide, 
refinement is. not difficult to make by introducing additional terms 
according to the 'character of the function. For instance, when circu- 
lation is present, the hydrodynamic functions for compressible flow 
transformed from those of incompressible flow for a sy mme tric body will 
give a negatively cambered airfoil. The necessity of applying correc- 
tion terms is obvious. 

This work was done at the Graduate School of Aeronautical Engineering 
of Cornell University under the sponsorship and with the financial assist- 
ance of the National Advisory Committee for Aeronautics. 


FUNDAMENTAL EQUATIONS AND THEIR PARTICULAR SOLUTIONS 


A steady, potential, and isotropic motion of a perfect gas in a 
plane satisfies a system of equations in the hodograph plane (refer- 
ences 1 and 2) : 


q 


dtp _ Po ^ 

dq P 



dqp _ Po 
d6> P 


q 


cw 

dq 


> 


J 


( 1 ) 


Here cp(q,0) and ^(q,0) are, respectively, the velocity potential 
and stream function; q and 0, the magnitude and inclination of the 
velocity vector, respectively; p and M, the density of the fluid and 
the local Mach number, respectively; and p Q , the value of p at q = 0. 
(See appendix A for definitions of all symbols.) 


This fundamental system, as shown previously (reference 2), can be 
solved and, in the case of ty(q,0), the particular integrals are 


0 VF v (T >{sin} VS 



~ ( 2 ) 
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|L[t v F v (t) cos v9 ] 

^*|t v F v (t) sin vej 


and 

• 6 and J (l - t)^ — 

and, in the case of <p(q, 0 ), are 


(1 - T )'VF v (T)l v (r)|_^jve 
(1 -T)-V v F_ v (x) Lv ( T )^jve 


> (3) 


(4) 


> (5) 


(1 ' t) '^ J^[ tVf v (t) *v sin ve ] 

(1 - t)“P ^|-|t v F v (t)s v (t) cos v ^| 




( 7 ) 


All these results are well-known except equations ( 3 ) and ( 6 ) which can 
be found in appendix B. Here F v (t) and F_ v (t) stand, respectively, 

for the hypergeometric functions F(a v ,b v ; v + 1 ; t) and 

F(a v - v,b v -~ v; 1 - v; t) for nonintegral parameters v. In case v 
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is an integer n, the definition of the second integral is modified. 
The new function will include a logarithmic term (cf. reference 2) and 
is denoted by F_ n (-r). The functions | v (t) and | (t) are defined 

fcy 


v M t ) = 2t l0ge tV ^ F v ( t) 


v| „(t) = 2t a log t" V / 2 F (t) 
dr e -v 


( 8 ) 


and the constants a v and b are defined by 


ay + b v = v - 3 


a v b v = 


- | Pv(v '+ 1) 


3 = 


7 - 1 


(9) 


and the variable t is 


t = -L sL 

23 c 0 2 


( 10 ) 


where 7 is the ratio of specific heats of the gas and c Q , the speed 
of sound at' q = 0. 

Furthermore, F v (t) is an analytic function of both variables t 

7-1 

and v. For fixed t < it has simple poles at v = -n and, 

7 + 1 v 

accordingly, possesses an expansion of the form (references 3 and ,4): 


f v<t> = f(r)T V •t i ^ i T'r(T)*'F 1 [Tl 


m=l 


(id 
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where 



T(t) 


- uf 



- x ) 1 / 2 + (l - ri 2 x) l/2 ] ri 


r ( a m) r ( m + 1 - »m) 

r(a m - m)r;(l - b m )ml(m - 1): 


, - f z + 1 

1 ^ 7-1 


> ( 12 ) 


This expansion is valid for all values of v except negative integers. 
When v = -n the limiting process gives 


00 8 


F (t) = f(T)T- n (x) - — Tm T m_n ( T ) F ( T ) 

u ~ — i -n + m 


m=l 



log e T 



(13) 


where the prime indicates that the term m.=n is to be excluded. 


From these particular solutions one can construct a solution for 
any boundary, in the hodograph plane. Once the functions qp(q,0) and 
ijr(q,e) are known, the flow pattern in the physical plane is given by 
integrating 


hx = 
he 

hy ■_ 
he ~ 



e _ _2. sin e 
he p 




cos 




( 1*0 
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and 


dx 

dq 


* q ( C ° ? 


Po . a dlA 

V 1 - — sm 9 rr- 

oq P oq J 


& -if. 

bq q\ 


sin 6 + — cos 0 

oq P 


( 15 ) 


where x and y are the coordinates of a point in the physical plane. 


MAPPING OF AN INCOMPRESSIBLE FLOW 


The potential flow of an incompressible fluid past a symmetric 
Joukowski airfoil without circulation is determined by a complex 
funct ion Q W( z ) : 


Q W(z) = £ + e + 


(1 + e ) 2 

i + * 



(16) 


where 


z = x + iy 

£ = 5 + ±r\ 


Here all the physical quantities are made dimensionless by dividing, 
respectively, every length and velocity by the radius of the generating 
circle and the free-stream velocity. Then € is a purely real number. 
The complex velocity w for such a flow is 


£ 2 U + 1 + 2ej 

(I + D(£ + e ) 2 


This equation expresses the relationship between w and t. It 
can also be interpreted as a transformation function from the £ -plane 
to the w -plane . In other words, if a flow field in the physical plane 
is given, its image in the w-plane is immediately known by means of 
equation (17). Since the present problem is formulated in the w-plane 
instead of the z-plane, the flow field in the z-plane can be calculated 
from the inverse transformation £(w), which is given by solving 
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£ 3 + (1 + 2e)£ 2 


e(€ + 2)w 
1 - v 5 



( 18 ) 


The three solutions of this equation are: 


^i(w) = P(v) + Q(w) - - y — 


^(v) = -|[P(w) + Q(v)] - ±|%v) - Q(w)] - 


1 + 2e 


(19) 


^3(w) = -^(p(w) + Q(w)] + ~[p(w) - Q(w)j - — 


+ 26 


where 


* w) ■ -? < 


( 1 . e) 3 , gliifi-Li > 

1 - V 


., 3/2 / 1 , o ^ 3 / 2 X / 2 

13 6(1 + 26 ) w 


(1 - w) 


3/2 


1 _ (1 - 0 3 (1 + 0 w 

(1 + 26)3 


> ( 20 ) 


Q( w ) = -i 4 


(1 e ) 3 


9e(l + e + e 2 ) 


1 - w 


i 3 3 / 2 6(1 + 2 6 ) 3 / 2 w 1//2 


(1 - w)3/2 


1 (1 - + 0 v 


(1 + 26)3 


This indicates that to a neighborhood in the w-plane there corre- 
spond, in general, three distinct neighborhoods in the £-plane or the 
z-plane (see fig. l) . The exceptional points where two of the three 
solutions become equal are given by dw/d£ = 0. This yields three 

points £ = 0, ty = - ^ ^ e , and £ = 00 about which two of three 

2+6 
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solutions join. Furthermore, it remains to investigate how the axis of 
reals, that is, v = 0,' is mapped in the £ -plane. By separating into 
real and imaginary parts, the velocity component v in terms of | 
and r) is: 

. _ 2 gT )(| 2 + e£ + T) 2 ) [(2 + e ) | + 1 + 2^ 

+ 21 + 1)(£6 + 2 el + e 2 ) 2 


Then v = 0 corresponds to 



(2 + g)| + 1 + 2e = 0 



In other words, the axis of reals of the w-plane corresponds to: (a) The 

axis of reals of the £ -plane, (b) a straight line through one branch 
point parallel to the T}-axis, and .(c) a circle. Moreover, while u 
increases from zero to one along v = 0, £^(w) decreases from the for- 
ward stagnation point to negative infinity, £g(w) increases from zero 
to positive infinity, and {^(w) decreases to a negative number. Evi- 


dently, branches I and II join smoothly at the line 


l + 2e 
2+6 


and 


branches II and III meet along the circle 



Branch III then lies entirely inside the circle of radius e/2. Hence, 
the flow over a body with a sharp trailing edge, complicated as it is,- 
remains essentially the same as that over the body which is symmetric 
about two axes with two stagnation points. Finally, by substituting 
equations (19) in equation (l6), the corresponding branches of the com- 
plex potential are 0 W]_(w) , 0 W2(w) , and qW^w). 


When circulation is present, having the magnitude required by the 
condition of finite trailing-edge velocity at an angle of attack a, 
the complex potential is 


0 W(.z) = (£ + e)e" ia + ^ + 


e' i ' a + 2i(l + e) sin a log e (£ + e) (21) 


i + e 
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with z = £ + and the complex velocity is 


w = 


£ 2 [(£ + e)e -ia + (1 + 

( £ + 1) ( £ + e) 2 


^ ia ] 


( 22 ) 


The inverse transformation function £(w) in this case consists in the 
solution of 


(e" la - v )£ 3 + [ee ia + (l + € )e ia - (l + 2 e)w ]^ 2 - e (2 + e)w£ - e 2 w = 0 

(23) 

The solution of this equation can, of course, be given, but the 
expressions are so complex as to make any practical manipulation dif- 
ficult. To bring out the complicated character of the Reimann .surface 
it is sufficient to determine the branch points and branch lines. Now 
the branch points are £ = 0 , £ = 00 and 

(l + e)(l - e 2a1 )^ 2 + e [3 + 6 + (l + e)e 2a ^C + 2 e[e + (l + e)e 2a £] = 0 


of which two are inside the body; and the wind axis in the hodograph 
plane, that is, the axis in the direction of the wind, corresponds, even 
for small values 'of a, to curves in the £ -plane: 

a(l + e)(! + 1 + 2e)r^ + e[(2 + e)6 + 1 + 2 <?|t) 3 -+ a(l + e)^ 3 + 

6(1 + e) I 2 - (3 - 2e - 6 2 ) | - (1 + e )^]t| 2 + el(6 + e.)[(2 + e)t + 

1 + 2ej t) + a( 1 + e)^ 2 (| + i)(| > e ) 2 = 0 

This equation is of degree four in i) and five in |. It would be 
interesting to trace, this equation for given values of e and a. 

This was not done for the sole reason that the exact problem is too 
difficult, if not impossible, to handle. It should be remembered that 
the present procedure of constructing a solution of a compressible flow 
depends primarily on the knowledge of an incompressible -flow solution 
in the form of power series. Even in the case of zero circulation, as 
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can be seen from equations (20), the series expansion of the complex 
potential yields sets of coefficients which are, at least, double 
infinite series (cf. appendix C). One can easily see that the initial 
task alone is a formidable one, not to mention the later computation 
of the compressible-flow solution. This ‘amount of labor is particularly 
unjustifiable iri this problem because, even if one starts with an exact 
solution, he may still end up with a body very different from what he 
expected because of the distortion due to compressibility. For this 
reason deviation from the formal procedure is recommended wherever it 
is convenient. In this report, the case of small thickness and small 
angle of attack is studied first. 

Before proceeding further, it is mentioned that the Joukowski air- 
foil is characterized by having a cusp at the trailing edge. Although 
the velocity is finite there, the acceleration is infinite and, con- 
sequently,. the analyticity of the mapping function breaks down; namely, 
d w/dw = 0. In order to preserve this property in the case of compress- 
ible flow, the stream function i(q,6) must satisfy 

= 0 

dq 

** = 0 

be 

In the case of symmetric flow, the first is satisfied identically. The 
second serves to determine the speed at the trailing edge. 


APPLICATION TO THIN SYMMETRIC AIRFOILS 


Consider the case of a thin symmetric airfoil present in a uniform 
flow without circulation, namely, a = 0. Since the thickness ratio of 
a Joukowski airfoil is directly proportional to the parameter 6, for 
thin airfoils e must be small. Suppose that the airfoil is so thin 
that terms higher than e are negligible. Equation (l8) is then 
simplified to 

f + {1 + 2 ^_J¥L. = o 



(25) 
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of which the solutions are: 


£l(w) = (1 + 2e)(l + 


M w ) = ~h (1 + 2€)(l - 


1 


1 - €1 w 


1 - w 


i 


1 - €1 V 
1 - V 


(26) 


where 


1 - 2e 
1 + 2e' 


being, of the order of magnitude unity. 


By making such an approximation it is seen that the third solution 
disappears. This is explained by the fact that the third branch lies 

wholly inside a circle of area ^ which, according to the present 

approximation, is zero; therefore £g(v) = 0. The branch points are 

now located at ^=-— (l + 2e) and C — °°, connected by the branch 
line (fig. 2) 2| + 1 + 2e = 0. 

The problem is indeed simple. The question, however, is: What 
flow is represented or, in other words, how is the complex potential 
modified? By using the simplified transformation function, a straight- 
forward integration gives 


Q W( z ) = £ + e + log e ^1 + 


(27) 


The flow in the £-plane is then that produced by the superposition of 
a sink and a source at, respectively, £ = 0 and £ = -2e in a uniform 
flow. It is very easy to see that, if e is small, the vanishing of 
the stream function corresponds closely to a circle with radius 1 + e 
and center at £ = - e. Hence for small values of e the curve in the 
z -plane will not differ much from a Joukowski airfoil. By substituting 
^1 and ^2 from equations (26) in equation (27) the two branches 
of the complex potential are 
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0 Wi( V ) = -i 



p 

1 - €1 w 

1 + (1 + 2e ) \ 


N 

1 - w 


1 + 2g 4e /f 1 - 

2e L IT^7 /\ 1 - w j + 


3 W 2 (w) = -i 


1 - (1 + 2e) 


1 + 2e 
2 e 


log. 


JlZj5 

V 1 - V 

1 - -li_ /- ( L ^i 2 -) 

1 + 2 e / \ 1 - w / 


L (28) 


1/2 


+ 1 


By taking e = 0.08, the holograph i|r o (q,0) = 0 from equations (28) 
and that of a symmetric Joukowski airfoil are compared in figure 3. It 
is seen that, except in the neighborhood of the singularity w = 1, 
the agreement is quite close even for a thickness ratio of about 
10 percent. 


For the case where a weak circulation is present, by the same 
method of approximation, equation (23) simplifies to 



(29) 


Here a is assumed to be the same order of magnitude as e, and there- 
fore the product of e and a is also dropped. The solutions of this 
equation are 


£l(v) = - 


1+ 2e 


2~\ 


1- «i v 


1 + U _l +a i i 

k ) l-v / [^1 -V U . v) l/2 (l . 6l 2 v )l/2 


(30) 


6 2 (v) --^Ml 


1 - W 




+ ai 


1 - v 


(1 -w) 1,2 (l- 
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the corresponding complex potential is 

0 W(v) = (1 - ai)£ + e + 1 + 2 g g + al log e ^1 + + 2ai l0 Se (5 + 2e) 

(3D 

The body in the £-plane is again a circle, provided both e and a 
are small. A substitution of ^ and £2 from equations ( 30 ) in 
equation (31) gives, respectively, 0 W 1 (w) and 0 W 2 (w). 

It should be noted that, in the previous case when circulation was 
present, it altered radically the nature of the singularity of the flow. 
Now by comparing equations (30) with equations (26) one can see that 
circulation has introduced essentially no complication and hence the 
problem can be similarly handled as in the case where no circulation 
exists. 


EXPANSIONS OF 0 W(w) 

When q=_0 . - For thin airfoils, as seen from the previous section, 
the domain of interest has two branch points £ = (1 + 2 e) and £ = <*> 

the corresponding points in the w-plane are w = 1 and w = — at 

61 

each of which 0 W(w) becomes singular. As a result, two different 
expansions are required to represent each branch of 0 W(w). In the~ 
case of 0 Wi(w), the origin is a regular point; for o^2( v ) if is a 
logarithmic singularity. 

Consider 0 W 1 (w) , for the case, a 

stagnation point | £| ^ 0 and qW-^w) 

» 

/ \ t* 1 + 2e 

0 Wl(w) = + 6 + 

>1 


= 0 first: About the forward 

can be simplified in the form: 


e(l + 2 e) 
Cl 2 


0(e 2 ) 


(32) 


^2 “ " 


2€V 
1 , - W 


By means of the relation: 
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Q Wp(w) can be reduced to two elements (l - w) ^ 2 ^1 - ep 2 w) "^ 2 and 

(1 - w) 1 / 2 ^ - € l 2 v) 1 / 2 ‘•which possess the following expansions in the 
respective domains: 


(1 - ^ S n (i) w n 

(1 - w) 1 / 2 ^ - e x 2 w^ 2 =£3 pS^^w 11 
v ' 0 


I w | < 1 


(33) 


where 


Sn (1) = (-l) n E 

m=0 


n 


lS„ (i) - (-1)“ 




2( h-m) 


'2(n-m) 


L (34) 


and 


(1 - ») -V*(l . n 2„)-l/2 . £ Sn <0) U 2n “ n - 


n s n (°) u -» 

,i 


(1 - w) l/2 (l - ei 2 w) l/2 = f 1 S n (°)e 1 2n w n + 

0 - 


1 <l w l < “2 

e l 2 


L (35) 


00 



16 


MCA TN 2356 


where 


Sn (o) = (-D n £f‘lY '2 Vi 
0 \ m /\m + n / 

- <-> n t(H * Vi 

0 \m/\m + n/ 


(36) 


A simple reduction then yields 


^(w) = -A 0 - A n w n + o(e 2 ) |w| < 1 


(37) 


where 


An - ^(s„ (l) - A-l (l) ) ♦ ^l>![sn + l (i) - (l ♦ ■. 1 8 K« 1) * 

f iVi (1) ] ♦ li ir^[-v 2 <1) ♦ (* * < 2 ) Vi (I) - 


26i 2 S n (l) - S n (1 > + 6l 2 S n .! (l >| 


(38) 


-Here' "A?) is~zero, as is easily shown from equation (27); in the above 
expansion it would have been of the order of e 2 and hence is included 
in the symbol o(e 2 ); The velocity potential and stream function are, 

neglecting 0 (* 2 ), 


j^l _ ~ A o ~ A n^ n 


cos n& 


00 

3+i = YL A n“l n sin ne 


q < 1 


(39) 


/ 
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If the expansions (35) are introduced, in the annulus region one will 
have 


qW^w) = B_ 2 w" 2 + B^v" 1 + B 0 +if (b v w v + C V W V ) 1 < |v| < -i- 


where 


( bo) 


Bn = - 


(l + 2e) ( 3 + 8e) 


B 


(1 + 2e) 

2e 


8s 

2 


B n = - Lk-t 2 e ) 1 


3 -2 


8e 


1 + 2*7 (o) o (o)) (1 + 2e) 2 „ ( 0 ) 

g \ n+1 •" V ) ~ j~ “ — l S n+l V 

(1 ; t 2<)3 (fl»l (o) -l S n +2 (o) ) 

cv = 4^(sn (0) - ^ 2 s n+ i< 0 )) . lSn (o) 


p _ 1 + 2e(o (o) r 2 0 (oh (1 + 2 s ) 2 0 

C l/2 i V s ° ' 6 1 S 1 / 47 I s 

(3. + 2g) 3 f g (o) _ p 2 g (o)j 
8s V- 0 111/ 


(o) 


> (4l) 


v = n + 


* J 



18 


MCA TN 2356 


Separating into real and imaginary parts,’ and are, 

approximately. 


Q cp^ = B 2 q -2 cos 20 + B ^q -J - cos 0 + B q + ^ ^B v q v ~ C v q" v J sin V0 


o+L = B_ 2 q -2 sin 20 + B_ ± 


E ( B v1 V - Cvl- 1 ") 

0 v 

q" 1 sin 0 + ^ (B v q v + C v q _v ) 


? (42) 


cos v0 


for 1 < q < 


e l 2 


In the case of Q W 2 (w), as w = 0 is not a regular point, a Taylor 

expansion does not exist. It can be expanded, however, in the following 
form: 


o W 2 (w) = ^ A_ n w _n + Y2 + °( e2 ) M < 1 

0 2 


(43) 


where A n is defined in equation ( 38 ) while A_ 0 , A_-p and A_ 2 are 


-o 


A , = 


l -l 


3(1 + 2 e) 
4 e 

(1 + 2 e ) 2 


(1 + 2 c ) 3 

4e 


(44) 


—2 

On the other hand, in the region 1 < |w| < e the expansions are 


Q W 2 (w) = B -n w_n " 1 HL ( B v wV + C v w " V ) + °( e ^ 


(45) 


1 



NACA TN 2356 


19 


Here the constants B_2, B_^> B Q , B v , and. C v are defined in equa- 

tions (4l). It is clear that by circling once about w = 1 branch I 
goes smoothly over into branch II, as is easily seen from equations (40) 
and (45). Similarly, the velocity potential and stream function are, 
omitting o(e^), ' 


o q>2 = 2, A -n t l~ n cos nQ + y A n q n cos n G 


0 < q < 1 


L (46) 


3^2 = YL A n q n sin + y A n q n sin n0 
1 2 


and for 


1 < q < € -2 


0^2 


= YL B -n q_n cos + (B v q V - C v q" v ) sin v6 


> (47) 


o^2 = |Z B - n q " n sin n6 ~ ( ] 


B v q v + C v q~v^ cos 


When a ^ .0 ♦ - When a is finite but small in the case of a thin 
airfoil, equations (30) show that 


= ^ (0) + aC x (l) 
= C 2 (o) + a^ 2 (l) 


(46) 


where ^ and are the transformation functions for a = 0 

and and £ 2 ^ are correction terms due to circulation. There- 

fore, they are defined as follows: 



20 


NACA TN 2356 



?1 <°> . .1^4 ♦ 




-1S, (1) — i 


1 - w u^(iT^ 


-i5 


(i) 


x - v (1 . V ) 1/2 (i . 


(^9) 


Similarly, the complex potential can be written as 

Wi = Wn ^ ° ^ + a Wi ^ ^ 

o w l o w l + a o w l 

c w 2 = 0 W 2 (o) + a 0 w 2 (l) 


' (50) 


where 0 W-^°) and <^2'°^ are ( i e i"i rie( i in equations ( 28 ) and 


and 0 W 2 ^ are found to be 


-i q vM = e x (D - * i|Jtl l°8e (h (0) - 2«) - 

1, , (o) 

26 1 ' ’ ^Co) ({!<«) ♦ 2e ) 

-iw 2 (i) = SgO) - {g(<>) + i-|±S I 0ge + 2e) - 


1 *- (o) 

^ l°8e fe 


(1 + 2e)t 


( 1 ) 


E 2 (o)(£ 2 (o) + 2e) • 


^ (5D 


which are valid if -f 0 
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To expand 0 W(w)^, "besides the fundamental forms, expansions- (33) 
and (35), it is required also -to consider 

l°g e [d " 2 e)(l " w) 1 / 2 + (l + 2e)(l - e^w) 1 ^ 

and 

1 - v) l/2 ♦ (1 - 6 /„) l/2 ] 

For the first expression let 

F(w) = 


dw 


log. 


d 


2e)(l - v) 


1/2 


(1 + 2e 


)(l _ ei 2 w ) 


1/2 



By direct differentiation one obtains 

F(v) = ~^ (1 - W ) _1 ^ 2 (l - e x 2 w ) _1 / 2 

By making use of equations (33) and (35), a term-by-term integration 
gives 


log e [(1 - 2e)(l - w) l/2 + (1 + 2«)(l -€i 2 w) 1/2 ] 


_^^fn-r__^ + 2 

2 hr n e 


U) 


% 


and, aside from a constant, 


M < 1 (52a) 


loge {(1 - 2e )(l - v) l/2 + (l + 2 € )(l - ei 2 v ) l/2 ] = 


Hi 

2 


£ y (“ s n+l <0)w ” V * Sn (o) 'l 2n » V ) 


(52b) 



/ 
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A similar method yields 
l°g e 


[(1 - v) l/2 + (l - = 


. i £ S2±. V» + log, 2 
? n , e 


(i) 


and 


log e 


[(1- w) 1 / 2 + (l - e] 2 *) 1 ^] = 


M < 1 


(53a) 


| loge w + iE ^(s n+ i (o) 


€i 2n+2 W V - S n 


( 0)4 


1 < |v| < e-f 2 (53b) 


( 1 ) 


Furthermore , 




can be written in an expansible torm: 


--H 1 + « n 
26 


1 - w 


1 1 1 - € l v 

Collecting all the individual expansions, one finds 

-i jl^ 1 ) = 2«i + E A n 1 ^ vI1 + °( e2 ) M < 1 


( 5 * 0 . 


where 


A n (1) - -(l + S n (i>) + ^-^(s n < lJ - *i 2 Vi (1) ) ' 

1 + 1/, ' 2, (i)^ 1 If, s (D\ 

-JT - SV 1 ' 5 1 S »-l / ' IT nV, 1 S n ) 


^(Vl (l) - Bn 11 ’) 


2e 


( 55 ) 



NACA TN 2356 


23 


(l) 2 

Here is of the order of e and therefore neglected. The real 
and imaginary parts of the correction term of Q W(w) are, to the order 
of approximation considered, 


+ £>n (l V 


sin n 6 


q < 1 


o^ 1 ) = y A n ^^q n cos nd 
2 


(56) 


_p 

In the region 1 < |w| < } with the exception of a constant 


-i = - 1 - loge we 111 + 


? I 1 + iT 


y + C v ^ 1 ^w _v ^ + o(e 2 ) 


(57) 


where 
B ^ 


(1) 


- Is (°) c 2 1 + 2e /s S (°Ac 2 1 1 + e. (0) 

- j-n +1 e l g V n+1 ’ n / 1 ~ Wv ' Sn 

Vi ( °V - 4^(Vi (0, n 2 - s n (o) )«i]«i 2n 

- V°> - . ^(ol) . ^(o) + 


>(58) 


_L. S (°) +-UL 2 £/ r s .,( 0 ) 
4 ev n 2 e 


(Sn + 1 (0) - Sn M )'l 
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The corresponding functions 0 <p^( q ) and Q \|r^( q ) are 


o‘ p l (1) ■ - e > - (l + j)l' n sin ne 


£ (Ev U) q V - C v (l) q- V ) 


sin v6 + Constant 


> (59) 


0 *i (1) = io ® e 1 + £ (* + ^) q_n c ° s ne + 

■ OQ y t 

> ^B v ^^q v + V j cos v9 + Constant 

Similarly, for the second branch the expansion is‘ 

-1 0 W 2 (i> = C^ 1 ' - i *1 - i Inge » + i c n U) ^ + 


|w| <1 

(60) 


where 


(1) _ l+4e 


= -1:7- l0Se e " 2^ ' 0Se 


1 + 26 


+ 2 


C n (1) = -(l - S n <i>). - Li^l(s n <i) - 

+ S ^ (1) ) - + S “ (i> ) + 
i -|j Si ( s n (1) - s n-l U) ) e l 


L (61) 
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For the region 1 < Jv| < , it is 


~ i o w 2 (1) = l0 s e w - lo e e wertl + 


00 / v 

£ i 1 + - 


E(b v < 1 >v v + c v w- v ).0 U 2 ) 


(62) 


where a constant of order unity is again left out. The corresponding 
real and imaginary parts are: 


^2 (1) = Yl (« " + YL C n (l) ^ n sin 


n0 


q < 1 


> (63) 


= C o (l) - ^ 1 °6 e < l + JZ C n^ 1 ^ Q n cos n 6 
dz 2 


and, aside from a constant, in 1 < q < ^ , 


00 

°9 2 (1) = -^ Q + (n - 6) - 2Z. ^1 + ~jq _n sin n 6 - 
f~ ( B v ^ 1 ^q V - C v ^ 1 ^q" V ) sin vG ■ 


^ ( 64 ) 


o*2 


(l) 1 + 4 e 


log e q + Y2. ^ + ^) q ~ n cos n6 " 


C(. 


B„ (l) q v + C, 


(D q -v) 


cos vQ 
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CONSTRUCTION OF A SOLUTION FOR A SYMMETRIC AIRFOIL 

The complex potential 0 W for a flow .over a thin Joukowski airfoil 
has now been transformed and expanded in power series of w for both 
the case of zero circulation and of a weak circulation. Knowing the 
form of expansion of Q W, a solution relating to a flow of incompres- 
sible fluid about a certain body' can be obtained im mediately by 
replacing each term of the series by the proper particular integral as 
given in the section "Fundamental Equations and Their Particular 
Solutions." For instance, the complex potential in the case of zero 
circulation is 

00 

■ qW-^w) = - T~ A n w n | w | < 1 

0 

The complex potential Wp(v,r) relating to a compressible flow is, 
accordingly, 

Wi(w,t) = - A n F n (r) ( T )w 11 + 7^ V c m(fl T ^ m V eme ' 

M < 1 <©> 

of which the imaginary part gives the stream function, namely, 

t]_(q,6) = ]m[w 1 (w,T)J (66) • 

Here the definition of F n ( r )(T) may be either F n ( t )/f ii (t 1 ) (refer- 
ences 1 and 2) or F n ( T >/ T-^^t-l) (references 3 and ^) , where t = 
when q - 1. It follows from equation (66) that the stream function is 

00 00 

= A n q I1 F r1 (r) ( T ) sin n0 + 2- Vm(fl T l)\i sin m6 " 

2 2 ' 

V 2 |;[(Fi T i) n *n ne ] n=2 


q < 1 (67) 
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From equation (5) the velocity potential is 

00 

9l_( < l> 0 ) = Aq - (l - t)"P YZ A n1 IlF n^ r ^ T ^n( T ) cos n6 + 

2 

cos m0 - A 2 c 2 ^ [^iT 1 ) n ^ n l n cos nej^ 

q < 1 (63) 

Both series (equations (67) and (68)) converge only within the unit 
circle q = 1. In order to have a function to cover the whole field of 
flow, both series must be continued across the circle of convergence. 

There have been two alternative methods proposed to effect this con- 
tinuation. The first method (references 1 and 2) is to have the solu- 
tion corresponding to the. flow in the upper z-plane arranged in one 
w-plane . After the series, for example, equations (39), (42) , and (46), 
are modified, they are required to be continuous with continuous normal 
derivatives on the circle q = 1 (fig. 4). On account of compressi- 
bility, the modified solution will be discontinuous in normal derivative 
and hence W 1 _ and W-l can no longer represent the same function. To 

correct this, one set of coefficients is left free and an extra series 
is added to W-^. These two sets of coefficients can then be determined 

by the condition at q = 1 (reference 2). The second procedure (refer- 
ences 3 and 4) is to proceed from a Taylor expansion, equation (67), say. 
Expressing the function F n (r) in the form of equation ( 11) gives rise 
to. a convergent double series which, when the order of summation is 
changed, is transformed into a form being convergent even for values of q 
greater than unity. Analytical continuation is thus accomplished. In the 
following discussion, the latter procedure is adopted for simplicity. 

At this point, the convention for- numbering the branch of Q W(w) 
is altered. Instead of dividing the z-plane fore and aft, it will be 
-more convenient to divide the plane into upper and lower halves. Then 
Q Wi(w) indicates the upper half and 0 W 2 (w), the lower. The parts 

corresponding to the regions q < 1 fore, 1 q < ?1 , and 0 < q <1 
aft are, respectively, 0 W-l_, 0 W x , and Q W 1+ . 

By the second method, replacing F n ( t) in equation (65) by equa- 
tion (11) and, for evident reasons, adopting. the definition 

F n^ = F n^ T V^1 ( T l) ’ e 9 u ^tion (65) can be written as 

00 00 00 

W].(v,t) = — f ( t ) YZ A n (tv) n + Y~ A n (tw) n )■ — 5S_ tVfJt) 

0 fee fe n + m 

00 

Vn,(ri T l) m ^ Ime - A 2 = 2 


+ 
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where t = T/Tp. The first series can he summed and, by changing the 

- order of summation, the second series can be transformed into a single 
series. Thus 


Wi(v,t) = f(r ) ^(w) - c m(f F 'l T l) m im ( T ) eim0 v m " 1 0 w l dw + 

m=l J 0 

- (69) 


^^2 

where w = tw and * T F m ( t). The first term corresponds to 

that resulting from the asymptotic summation (references 1 and 2) and, 
in the special case 7 = -1, reduces to the Karman-Tsien approximation. 
The integral in the second term is bounded and the series can be shown 
to be convergent for q > 1 (references 3 and 4). 

By taking a path starting from the forward stagnation point to an 
arbitrary point in the annular region, namely. 


Ow p® 

J w m " 1 0 W 1 (w) dw = J w m_1 0 W 1 (w) dw + J w m-1 0 W 1 (w) 


dw 


where ^ is a fixed point in the annular region, a simple transforma- 
tion leads to 


W-l = JY B -n F -n (r) ( T ) w_n + 1 Y B v F v ( r )(x)wV + C V F_ V ( r > (r )w- v_ | 

0 0 L — 1 


( -A 2 c 2 + B -2 p 2)^[(^l T l) n ^n ein ^] n _ 2 


(70) 
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by replacing Q \l± by the proper expansion under the integral sign, 
where 



Vi v+v 


CyV V+W \ 

-v + m / 


b -2 lo Se W 1 


(71) 


The integral in equation (71) 


can be evaluated in the £ -plane; namely, 



dw 





-1 w . dW 
01 d£ 




The fact that 

on 

Equation (70), evidently, is the continuation of equation ( 65 ). Pro- 
ceeding further to a fixed point w 2 inside the unit circle but below 

the real axis, by describing a loop surrounding the point w = 1, the 
following proper transformation is arrived at: 


[(m) n v lne ] 


is a solution is shown in appendix B* 


(r), 


*1 = fl Vm(fi V*"* * IZ A_2F_ n Uj (T) w“ n + 
2 ' 1 


£ A n F n (r, < 


t)w 11 


(72) 


where 


»m(?2) * 


<2 




2-- An»2' n ' M ° f- KPs** 

-n + m rT n + m 


A _2 lQ g e v 2 


(73) 
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The path is quite arbitrary; in order to preserve the property of sym- 
metry of the function it is so chosen that the s are real. Equa- 
tions (65), (70), and (72) together represent the first branch of the 
Riemann surface covering the entire upper z-plane. Separating into 
real and imaginary parts, the stream function is 


tl(q.e) = ^|B v q v F v ( r )(T) + C y q- V F_ V ( r ) ( t)1 cos V0 

0 L - J . 

^B_ n q" n F_ n ( r )(T) sin n9 - ^AgCg - 


"«] 


1 < q < er 2 (74) 


and 


00 ^ , 

^l+Cq^) = h m(^2) c m(f r i T l) m %i( T ) s in m9 - 
2 


^ A n q n F n ^ ( t) sin n0 + 
2 


f;.A_ n q-Vj^(T) sin ne 


q < 1 


(75) 
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The velocity potential corresponding to equations (7*0 and (75) is 


<Pl(q> 0 ) = (1 t) 


-3 


B_ n q" n F_ n ( r )|_ n cos n0 - 


T_ [B v q v F v ( r )(T)l v + C v r V F. v (r) ( T )l. v ] Sin V0 + 


(^2 C 2 - B - 2C2 )“;((fl T l) n *^n cos n «] n=; 


L 1 < q < € X ' 2 (76) 


and 


Ti + (q, 0 ) = (1 *■ t)-^ 


- IlI (jpi T l) m ^m( f ) £ni cos m0 - 


Y~ A_ n q _n F_ n ^ ( t ) |_ n cos n0 + 

r 


^Z'A n q n Fj r )( T )f jn 


cos n0 


0 < q < 1 (77) 
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By taking a path from the forward stagnation point hut in the 
reverse direction, the second branch of W(w,t) can be' joined. Nov 
the solution 0 W(w) to be ‘considered in the integral- appearing in 
equation ( 69 ) will be 


|w| < 1 


3 w 2 _ - - v 
0 


2 CO . \ 

d W 2 = H B -n w_n - 1 TL ( B v wV + c v w " V ) 1 < M < e l" 2 
0 O' 


> (78) 


2 

o W 2+ = H A -n w_n + I_ V +n 
1 0 

Complete symmetry with respect to the u-axis makes it 
repeat the process. 


|w| <1 - 


J 

unnecessary to 


In the case of the flow with circulation, the complex potential 
consists of two parts as shown in equations (50). The one which corre- 
sponds to zero circulation is given in equations (37), (4o), and (43). 

The functions for the compressible flow are then given from equations ( 67 ) 
to (77). The part which takes into account the effect of a weak circu- 
lation is 


- 1 o«l- (1) 

0 


I w | < 1 


-i oW^ 1 ) = B q - j^ .we* 1 * t (i + ^j v ~ n + t ( B v (l)wV 


( 1 ) 


w- v )- 


^ (79) 


log e we 


Ttl 


00 

- 1 o W l+ 1 ■ "27 log e w + XI c n l)wI1 

0 


0 < J w | < 1 
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By the convention adopted, > and ’ 0 w i+^ form' the 

first branch of and correspond to the whole upper z-plane. The 

same method yields a solution for a compressible flow: 


-iw 1 _ (l) = Y A n (l) F n ( r )(T)w n + Y~ $*K eim9 ~ 

0 2 . 


-iW 


l (l) = Bo - T- U “ 8)1 - T- r 


T 


' £ (l ♦ tyJ'V* ♦ £ £<», W.V + 

C v (l) F. v (r)w "^] (« ~ e )i + J (1 - t) P ^ 


f(80) 


CO 

-iWl + (l) = Yi V ( 2 i %(^lTl) m t m (T)e me - 


f l 1 + s) C » |r[('Pl T l) n +n elne ] - 






L c n (l)p n (r| «" 


0 
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The real and imaginary parts are 


+ 1 - 


.An^^n^^T) cos n© + 


tl 1 )c m('P F l T l) m +m 


cos m0 - 


n 

YZ ( X + i) c n ^[(f F l ,r l) n in cos n0 ] 


q < i 


*l (l) - - Hr r (1 - T)e T- - £(* + |)^-n (r) ( T > 

YZ (Rv (l) q V F v (r) (T) + C v (l) q- v F_ v (r) (x[l cos V0 + 


cos n© + 


>(81) 


Constant 


1 < q < e- 


♦u (1) - - i 


PT 


2 € 


(1 - t) P + YZ h m (2)c m(V T i T ]) m V cos m0 ‘ 


£- i 1 + n) C " s[(fl T l) n *n CO. ne] 


£ C n ( 1 >qX( r >(T) 


cos n© 


0 < q < 1 



NACA TN 2356 


35 


CP 1 _ ^ - 1 ^ = (1 - t)"P < / Sin n 0 + 


tl 


sin m0 - 


_ 

ZZ ( 1 + n) c n jfo {('r r i T l) ni < f n^n sin n0 | ► 


= (1 - t)"* 3 j^r (l + |) < l” nF _n^ r ^ T )^-n 6in n0 + 

« 

00 

y~ [B n ^ 1 ^q V F v ( r )(T)| v + C v ( 1 )q- v F_ v ( r )(T)5_ y J sin V0 


Constant - ■— -' !' (« - 0) 

4 e 


1 < q < e- 


, (1) = (1 - t) 


00 

' P - Y^- h m (2)c m(V T l T l) m %i sin mQ - 


Y_ (l + |)<=n ^[(flTl^tnSn sl “ ^ + 

Y~ C n ^ 1 )q n F n ( r ) ( T )t n sin n0 - (it - 0) 

0 2€ J 


q< 1 
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In a similar fashion, the second branch can also be connected by- 
taking a path in the reverse direction and introducing successively the 
functions: 


-1 0 W 2 - (1 > = 


|v| < 1 


-i 0 W 2 


(l) 1 + - ni 

£7- log e ve + 


£ l 1 + s) w ' n 


) (b^/ 1 ^ + C v ^ 1 ^w _v ) 


+ Constant 1 < |w| < 


( 83 ) 


-i Q W 2+ (1) = -±- logg w + £: c n (l) w n | W I < 1 

.. 2C . . . ° J 

However, it should he added that, since the functions and W^ 1 ) 

do not merely differ by a sign, the h m ^^'s in W 2 +^^ and W-]^ 1 ) 

will not generally be equal. As a result, the two branches may not be 
joined with each other across the branch line y = 0, x>> < x < oo. In 

order to avoid this difficulty, the lower limit in the integral, 'for 
instance, equation ('73) > should be adjusted to make them equal. 


Having the stream function and the velocity potential so chosen, 
the flow pattern in the physical plane can be calculated by integrating 
equations (l4) and (15) to give the coordinate functions x(q,0) and 
y( v q,0) : . When the flow is ;free of circulation) the stream function is 
antisymmetrical about ‘the x-axis. By taking x = x^, y = 0 . as the 

stagnation point and x = Xr>, y = 0 as the trailing edge, the constants 

of integration can be determined on the unit circle for Q > 0. By 
symmetry, the reversal of the sign of 6 gives the lower half of the 
flow. 


When circulation is present, the symmetry property is destroyed 
and hence both branches have to be calculated separately. Since 
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both w( ^ and are properly joined functions, the determination 

of integration constants on the unit circle for both branches would be 
sure to give a pair of continuous coordinate functions. ' 


Cornell University 

Ithaca, N. Y. , August 10, 19^+8 
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<p 
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Co 
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APPENDIX A 
SYMBOLS 

magnitude of velocity vector 
inclination of velocity vector 
value of p at q = 0 
density of fluid 
local Mach number 
x stream function 
velocity potential 
hypergeometric functions 
nonintegral parameter 
integers 


parameters of hypergeometric functions 

functions defined by equations (8) 

ratio of specific heats of gas 
speed of sound at q = 0 
defined by equations (12) 
defined by equations (12) 
defined by equations (12) 
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*1 


i 


y + 1 

7-1 


x,y 

0 W(z) 

z = x + iy 
€ 
w 

P(w),Q(w) 

u,v 

5(v) 

Ci(w),^2( w )^3( w ) 


0 Wi(w) , Q W 2 {v) , Q U^v) 
a 

* _ 1 - 2e - 

fc l 

x 1 + 2e 

0 W(w) 

*0 

<Po 

s n< ±) ^l s n { 1) 

An 

0^1 >0^1 


coordinates of point in physical plane 
complex potential 

geometric parameter of body 

complex velocity defined by equation (17) 

defined by equations (20) 

velocity components 

transformation function defined by 
equation (l8) 

solutions of equation (l8) given by equa- 
tions ( 19 ) 

branches of complex potential 
angle of attack 

complex potential for incompressible flow 
in w 

stream function for incompressible flow 
velocity potential for incompressible flow 

defined by equations ( 3*0 

defined by equations ( 36 ) 

coefficient defined by equation ( 38 ) 
defined by equations (39) 



4o 


NACA TN 2356 


t (o) f (o) 

*»1 '*>2 

t (1) t (1) 

y. (°) y-(°) 

o"l >o"2 

y (l) y 0 U) 
o"l >o"2 


B (1) C (1) 
D \) 


o'Pl (l) 'c,*l <1) 

Cb< l WD 

o^W 1 ’ 

W 1 (w,r) 


defined, by equations (4l) 
defined by equations (46) and (4-7) 

transformation functions 

correction terms due to circulation 

defined by equations ( 28 ) 

defined, by equations (51)- 

defined by equations ( 58 ) 

defined by equations (59) 

defined by equations (6l) 

defined by equations (63) and (64) 
defined by equation ( 65 ) 


F„< r >(r) - »„<t)|-Ti(t]) 




value 

of T 

at 

q ■ 

1 

value 

of T 

at 

q = 

1 


branches of stream function 
branches of velocity potential 


^1- *^1*^1+ 
t = t/t x 
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qWi(w) 


defined by equation (69) 


$ - tw 

* m ( T ) = 


Wi (l), y, (l), y, (l) 
o”l- >o"l >o w l+ 

o w 2- (l) .o W 2 (l) >oW 2+ (l) 


defined by equation (71) 
correction terms due to circulation 

defined by equations (8l) 

defined by equations (82) 


defined by equations (83) 
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APPENDIX B 


PABTICULAR SOLUTIONS OF FUNDAMENTAL EQUATIONS OF MOTION 


If t m e ine is a particular integral of differential equations (l), 
-^-(^ n e in0 ) is also a solution. For the differential equation then 


becomes 


ffc q . n2 ^<1 . „ 2 )* n = 0 

dq^p -dq ) P ' 


d_|^o q d 
dqL P dq\dn /J 


- St'fr - m2 )^ = ^ - m2 K 


As T|r n ( q) satisfies the first equation, the second equation determines 
n. Now 




By substituting ^jt n (<l)e in6 [] for it is easy to show 


dcp 


= 1 a/_- M + 0 fd + i M ^ e ine 

n Sqy P dndq P Sq n P dq / 

1 d_f_. foi ^ n + e ££ + i £ 0 * ^V ine 

n V* P SnSq P Sq n P dq / 


dq + 


Therefore 


9 = -i 


P 0 <1 


A(i i*B. eA 

chi\ n dq / 
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APPENDIX C 


SERIES EXPANSION OF COMPLEX POTENTIAL 


The Taylor expansion of the complex potential of a flow over a 
symmetric Joukowski airfoil is, for 0 < |w| < 1, 


where 


- C n + <X + 0 Dd 


, = C v + (1 + 0 D v 


00 

o W l = A n w11 

o w 2~|£ A “^-?£v V 

.W 3 = .if A n vH + f f B V .V 




n -J, 2J .k E 2 )^! - + «)» 


(1 + 2e) 


3{j-k+D 


v £ £ £®U ; X* ; % h 

+ 2e) \ 2 ' 


(1 + 2€) 


2(1 


^ £ t £ @m. 1 — 
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Figure 3.- .Zero streamline of incompressible flow around symmetric airfoil in hodograph plane. 

g = 0.08. 





